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FAST QUANTUM MECHANICAL INITIAL STATE 
APPROXIMATION 

GOVERNMENT INTERESTS 

This application discloses an invention made with government support under 

Contract No. F30602-01-0523 awarded by US Air Force, Air Force Material Command 

Air Force Research Laboratory/IFKF. The government may have certain rights in the 

invention. 

FIELD OF THE INVENTION 
This invention relates to quantum computing and to methods and systems to 
efficienfly calculate eigenvalues and eigenvectors Hermitian operators and quantum 
mechanical evolution operators with quantum computers and methods and systems to 
efficientiy calculate an approximate quantum state to be used as an input in a quantum 
mechanical system. 

BACKGROUND 
Intuitively, quantum mechanical problems offer great potential for quantum 
computers to achieve large speedups over classical machines. An important problem of 
this kind is the approximation of an eigenvalue of a quantum mechanical operator. In a 
recent paper published in 1999 in Physical Review Letters (Vol 83, p. 5162) and hereby 
incorporated by reference, Abrams and Lloyd present a quantum method for doing this. 
Their method is exponentially faster than the best classical method, but requires a good 
approximation of the corresponding eigenvector as an input. 
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There is currently a continuing need for a method and system for efiBciently 
computing a good approximation of the eigenvector as an input to the Abrams and Lloyd . 
quantum method. 

There is also a continuing need for a method and system for efficiently computing 
a good approximation of a quantum state (not limited to eigenvectors) as an input to a 
quantum mechanical computer or computation. For example, one would like to compute 
an approximate input to the quantum simulation algorithm. The quantum simulation 
algorithm is described in the book Quantum Computation and Quantum Information, by 
M. A. Nielsen and L L. Chuang, Cambridge University Press, Cambridge UK (2000). 

SUMMARY OF THE INVENTION 

The present invention is a system and method for use on a quantum computer to 
efficiently prepare the initial quantum state required by Abrams and Lloyd's eigenvalue 
approximation method. The system and method of the present invention is used to 
prepare a quantum register with an approximation of the eigenvector that is guaranteed to 
be sufficiently good to be used as input to the Abrams and Lloyd method. The present 
invention can be used when solving continuous Hermitian eigenproblems, e.g. the 
Schrodinger equation* on a discrete grid. 

Beginning with an eigenvector for a problem discretized on a coarse grid, the 
system of the present invention efficiently constructs, quantum mechanically, an 
approximation of the same eigenvector on a finer grid. This eigenvector approximation is 
suitable as the initial state for the eigenvalue estimation method of Abrams and Lloyd. 

Similarly beginning with a vector (i.e., a quantum state) for a continuous problem 
discretized on a coarse grid, the system of the present invention efBciently constracts, 
quantum mechanically, a vector (i.e., a state), which is an approximation to the 
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corresponding vector on a finer grid. Our system efficiently extends a vector of low 
dimension to one of high dimension, which is then presented as input to some quantum 
computation method, e.g., the quantum simulation algorithm. 

The features and advantages of the present invention will be more readily 

apparent and understood fix>m the following detailed description of the invention, which 

I 

should be understood in conjunction with the accompanying drawings appended to the 
end pf the detailed description. 

BRIEF DESCRIPTION OF THE DRAWINGS 

Figure 1 is a chart illustrating the steps performed on the various quantum 
registers according to an embodiment of the present invention. 

Figure 2 is a chart illustrating the steps performed on the various quantum 
registers according to another embodiment of the invention. 

DETAILED DESCRIPTION 

For purposes of illustration only, and not to limit the scope of the present 
invention, the invention will be explained with reference to the embodiments of the 
invention indicated in the drawings. One skilled in the art would understand that the 
present invention is not limited to the specific examples disclosed and can be more 
generally applied to other initial state preparation methods and systems than those 
disclosed. 

The key component in the Abrams and Lloyd method is quantum phase 

f 

estimation, which is a method for approximating an eigenvalue of a unitary matrix. 
Quantum phase estimation is also described in the above referenced book of Nielsen and 
Chuang. We give a brief outline of this method below. 
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Let Q denote a 2*" x 2*" unitary matrix. We want to approximate a specific 
eigenvalue of Q. Phase estimation does this using the corresponding eigenvector as 
input. The Abrams and Lloyd method deals with the case when this eigenvector is not 
known exactly. Referring to Figure 1, consider a quantum computer consisting of three 
registers 140, ISO, and 160 with a total of b + m-^w qubits. The first b qubits in register 
150 are all initially in the state |0>. The second register 140 with m qubits is initialized to 
some state | which must approximate the eigenvector in question sufficiently well, as 
will be seen. The last w qubits in register 160 are work qubits for temporary storage. 
The w qubits are not important in our discussion here, and we generally omit discussion 
of them below. 

Since Q is unitary and therefore normal, the state | tf^ can be expanded with 
respect to eigenvectors of Q. Omitting discussion of the work qubits in register 160, the 
initial state of the algorithm is 




where |«> are the eigenvectors of Q. Placing the first register 150 in an equal 



superposition, using b Hadamard gates in step 170, transforms this state into 




Next, powers of Q are applied in step 170 to create the state 
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26-1 

Since is unitary, its eigenvalues can be written as where (pu €R. We can 

assume that cpu e [0 J) and consider the approximation of one of these phases instead of 
the approximation of one of the eigenvalues. Equation (3) is equal to 

r ^tm^ 

It is easily seen that the inverse Fourier transform performed in step 170 on the first 
register 150 creates the state 



( i: K¥'u.i)i?> j i«>, (5) 



(6) 



where 

In step 180, a measurement of the first register 150 produces outcome j 190 with 
probability 

and the second register 140 will collapse to the state 

y;£«£fa|)|,,), (8) 
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represented by the register 200; register 210 contains the work qubits after the 
measurement 180 as known to one skilled in the art. 

We remark that for special case when the eigenvalues cpu can be represented 
exactly with fr-bits (i.e., 2^9 u is an integer), equation (5) simplifies to 

When the eigenvalues are of this form and are distinct, a measurement in step 180 of the 
first register 150 will cause the second register 140 to collapse exactly onto the 
corresponding eigenvector in register 200. 

Recall that the system and method of the present invention are to achieve an 
approximation of the phase that corresponds to an eigenvector \u') using a quantum 
computer, that the state |^) is an approximation of this eigenvector, and that the 
eigenvalue is obtained from the value of the outcome j 190 by e^^*-''^''^ is of the form and 
approximates e^"* For instance, one is often interested in the eigenvalue 
corresponding to the ground state or in low order eigenvalues- We define A(9o, 9i) = 
min ^ |jc + - ^J'ol }» 9^o.^i ^ ^ ^® firactional part of the distance between <po and 
cpi). Then a measurement of the first register produces an outcome from the set G = 
{/:A(j72*,^„.)<ik/2*,fc>l } with probabUity 
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and when = 1 the probability to obtain j such that a(j72*,^7„.)^ 2""* is bounded from 

below by — yl^a p • I n^^st be chosen in a way that this probability is large or 

preferably greater than Vz, which implies that has to be sufficiently large. For one 
embodiment of the present invention to obtain an approximation of q>u* with accuracy 2 " 
and probability at least \djl^ (l — e)^ equation (10) shows that the number of qubits h in 
the first register 150 must be 

l, = n+[log(n-^)]. (11) 

Quantum phase estimation can be used as an efficient subroutine to find 
eigenvalues. Consider a Hermitian operator H. The operator G(t) = e'^^ is unitary and 
has the same eigenvectors as H. We assume that G can be implemented efficiently and, 
therefore, can be used as the unitary operator in the phase estimation algorithm. For 
example, when H is local, i.e., it can be written in the form , where each Hj acts 

only on a small number of qubits, then G can be implemented efficiently. However, 
locality is not a necessary condition for efficient implementation. Lideed, G can be 
efficiently implemented for a many-particle quantum mechanical system with a non-local 
H. One skilled in the art will understand that it is possible to implement G for a wide 
class of non-local Hamiltonians. 

The Hermitian eigenproblem described above is solved on a discrete grid. One 
embodiment of the present invention addresses the case in which the grid is extremely 
fine. Clearly, a fine grid requires a large vector for the representation of the initial state 
of the algorithm. In general, it may not be possible to efficiently prepare an arbitrary 
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quantum state in a space with a large number of qubits. However, the present invention 
includes a method for the efficient preparation of an initial state. 

In one embodiment of the invention, the operator possesses an eigenvector for a 
coarse grid discretization of the problem, which was most likely obtained classically 
since the size of the problem is small, although one skilled in the art will understand an 
eigenvector obtained by any coarse method can be employed without diverging from the 
scope of the invention. Using this eigenvector, we efficiently construct an approximation 
to the corresponding eigenvector for a fine grid discretization of the problem. We use 
this approximation as the initial state of the eigenvalue approximation algorithm. We 
describe our method for a one-dimensional continuous problem on the interval [0,1]. 

Let if be a positive Hermitian operator, defined on a Hilbert space of smooth 
functions on [0,1]. Let VkirX ^ = 1»2,- . denote the eigenfunctions of H, ordered 
according to the magnimde of the corresponding eigenvalues; and without loss of 
generality we assume that 

jf' - 1. (12) 

Suppose that Hm is a discretization of H with grid size = 1/(1 + N). Let 

I J/ft ), /: = 0, 1 iV - 1 , denote the normalized eigenvectors of Hn , ordered according 

to the magnitude of the corresponding eigenvalues. The expansion of the Jk-th 
eigenvector in the computational basis can be written as 

8 
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Let I = 53 1^ iU + l)^jv )\j ) be the sampled version of Vjt( ) at the discretization 
points. Consider problems such that the eigenvector of interest satisfies 

1 



- Oihj,), (14) 



where ^ > 0 is the order of convergence and )||^ = ^ \xj\^ for 

I X) = 2iIo ^ -^j]-^ )' example, these conditions are satisfied when, for example, 
we are dealing with second order elliptic operators. 

Now, assume that the eigenvector [c/j^^"^) of H^^ has been obtained classically. 
This vector is placed in a log No qubit register 110 (see Figure 1). For N = 2Wo, we 
construct an approximation |f7j[^^) of |c/jt^^) by appending s qubits in register 120, each 

qubit in the state |0) , to Iz/i^**^) and then performing in step 130 a Hadamard 
transformation on each one of these s qubits in register 120, i.e. 



where fU) = ljf2'}. The effect of/is to replicate the coordinates of |c/;[^o^)2* times. 
According to the present invention, is used as input to the eigenvalue and 

eigenvector approximation method. When the result of the method is measured 

9 
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will collapse onto a superpositioii of eigenvectors according to equation (8). The 
magnitude of the coefficient of \u^^^^ in this superposition can be made arbitrarily close 
to one by appropriately choosing Nq. 

Consider two different expansions of |t7][^^^ : 

«>r> = ^<&io{'^). (17) 



The first expansion is in the computational basis and the second is with respect to the 

I P 

eigenvectors Hi^. We call d^^^ the probability of success. Equation (17) can be 
rewritten as 

f^fift <18) 

Taking norms on both sides and using (13) and (16) gives the inequality 

= l-lO'- (19) 

We will now bound (19) from above, and thus the probability of failure. The 
definition of \ui'*^) implies 
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2 W-l 



ic^n-it^r>)ir=E 



where Zj=;K]r=o(/4')and S^K/Oof =2^0^^) by (14). Applying the triangle 



2 

, (20) 



+ 0(AS,,). (21) 



inequality, we get 

«ifc((i-H)ftw) Vifc((/(i)-H)fejv<>) i \ 

SIli^^'^^Mi v^iiiv^-'>>ii I ; 

The definition of \vi^^) and the fact that ||vj^||^ = 0(l) imply that 
|v;^^> )[|=V7i/(l + 0(&^)). Hence, the sum above is equal to 
1 

E MU + mj^Kl + 0(fcivr)) -vU(f<J) + l)fc/yb)(l + 0{AiVb))P- (22) 
Since Vfc(-) is continuous with a boimded first derivative, we have that 

where x^ j -{j'¥l)h^ and Xj^^. =(/(7) + l)/i^^,jr = 0,...,iV-l • Clearly 

1^2,; "^ijI - ^{Ko )• Using (22), (23) and the triangle inequality, we obtain from (21) 

that 

Hence, the probability of failure is bounded from above by o{Nq^^^'^^^). It 
depends only on the order of convergence to the continuous problem and the number of 
points in the classically solved small problem. We can select an No such that the 
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probability of failure is less than ¥2^ no matter how much larger N is. By choosing a large 

we can make the discretization error arbitrarily small. Equation (24) implies that the 
probability of obtaining the eigenvalue e^^'^* with accuracy 2'^ is at least 

We remark that any classical numerical algorithm that computes an eigenvalue, 
satisfying a specific (nontrivial) property, of slNxN unitary matrix takes time Q(iV)- For 
example, one may want to find the eigenvalue that corresponds to the ground state. This 
j is true even if a matrix is sparse and regardless of whether the algorithm is deterministic 

' or randomized. It is merely a consequence of the fact that the algorithm needs to 

consider all the (nonzero) elements of the matrix, and there are at least Cl{N) such 
1 elements. Alternatively, in the restricted case when the matrix is diagonal finding one of 

I its elements is a problem at least as hard as searching an unordered list. The lower bound 

j for searching yields the lower bound in our case. 

I 

I In conclusion, our method provides a highly efficient preparation of initial states 

for eigenvalue approximation, requiring only a small number of Hadamard gates. Thus 
the method of Abrams and Lloyd, using the initial state prepared by the system and 
method of the present invention, computes the eigenvalue exponentially faster than any 
classical algorithm. The method of the invention can be generalized to higher 
dimensional continuous problems. 

In another embodiment of the invention, if we possess a vector that corresponds to 
a coarse discetization of a continuous problem then, under suitable conditions, we can 
efficiently extend it to a vector that approximates the corresponding vector (i.e., a state) 
of a fine discretization. Referring to Figure 2, we first place the original or given vector in 
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register 310* Assuming that the vector has dimension Nq this register has log No qubits. 
For a iV = 2Wo, we append to register 310 s qubits, in the state |o) , in register 320. Then 

in step 330 we apply the Hadamard transform to the appended qubits. See Equation (IS) 
and the explanation of the effect of the replicating function/ register 340 we have the 
combination of the two registers 310, 320, register 340 containing the approximation 
corresponding to a vector (i.e., state) of dimension N = 2Wo. This requires log N = log A^o 
+s qubits for its quantum mechanical representation. Step 350 represents a quantum 
mechanical system using the approximation obtained in register 340. Step 360 represents 
the final state of the system 350. 

Having described the embodiments of the invention, it should be apparent that 
various combinations of embodiments may be made or modifications added thereto as is 
known to those skilled in the art without departing from the spirit and scope of the 
invention. 
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FAST QUANTUM MECHANICAL INITIAL STATE 



ABSTRACT 

A system and method efficiently prepare the initial state of a quantum computer 
required by the eigenvalue approximation method of Abrams and Lloyd. The system and 
method can be applied when solving continuous Hermitian eigenproblems, e.g. the 
Schrodinger equation, on a discrete grid, and allows for efficient calculation of their 
eigenvalues with quantum computers. A system and method efficiently prepare an 
approximate initial state (not limited to eigenvectors) of a quantum computer required by 
a quantum algorithm as input. 
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